ntricity or obliquity tides have been proposed as the missing energy source that may explain 
the anomalously large radius of some transiting "hot Jupiters". To maintain a non-zero and large 
obliquity, it was argued that the planets can be locked in a Cassini state, i.e. a resonance between 
spin and orbital precessions. We compute the tidal heating within "inflated" close-in giant planets 
with a non-zero eccentricity or obliquity. We further inspect whether the spin of a "hot Jupiter" 
could have been trapped and maintained in a Cassini state during its early despinning and migra- 
tion. We estimate the capture probability in a spin-orbit resonance between ~ 0.5 AU (a distance 
where tidal effects become significant) and 0.05 AU for a wide range of secular orbital frequen- 
cies and amplitudes of gravitational perturbations. Numerical simulations of the spin evolution are 
\Q > performed to explore the influence of tidal despinning and migration processes on the resonance 

■ stability. We find that tidal heating within a non-synchronous giant planet is about twice larger than 

■ previous estimates based on the hypothesis of synchronization. Chances of capture in a spin-orbit 
resonance are very good around 0.5 AU but they decrease dramatically with the semi-major axis. 
Furthermore, even if captured, both tidal despinning and migration processes cause the tidal torque 

_-, ' to become large enough that the obliquity ultimately leaves the resonance and switches to near 0°. 

Locking a "hot Jupiter" in an isolated spin-orbit resonance is unlikely at 0.05 AU but could be 
J> ' possible at larger distances. Another mechanism is then required to maintain a large obliquity and 
create internal heating through obliquity tides. 
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1. Introduction 

During the past decade, one of the most surprising findings was the discovery of several extrasolar 
planets with periods down to 1 day and orbiting their parent stars at a distance lower than 0. 1 AU 
(e.g. Santos et al. 2005). Up to now, fourteen of them have been detected by the transit method. 
The coupling between radial velocity and photometric measurements allows direct and accurate 
determination of both their masses and radii. However, the comparison with theoretical evolution- 
ary models of Jovian-mass planets revealed unexpected conflicts between measured and predicted 
radii. The first to be discovered, HD 209458b, was found to be ~ 20% larger than expected (e.g. 
Bodenheimer et al. 2001; Guillot & Showman 2002; Baraffe et al. 2003). Since then, new compar- 
isons with models that include the effect of strong irradiation from the parent star (Baraffe et al. 
2003; Chabrier et al. 2004) indicate that three other transiting "hot Jupiters" may still display an 
anomalously inflated size and low density compared to Jupiter. Their characteristics are summa- 
rized in Table 1 . 

Many scenarios that have been proposed to explain this discrepancy invoke a missing energy 
source that would slow down the gravitational contraction and cooling of the planet. Guillot & 
Showman (2002) suggest that a conversion of only 1% of the stellar radiative flux into thermal 
energy at depth would be sufficient. The contribution of tidal heating due to a planet's eccentric 
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Table 1. Characteristics of inflated hot Jupiters discovered so far and the rate of obliquity tides 
heating within them for various values of the obliquity e. 





Rate of tidal heating 


(erg.s ') 


Name 


*Age[Ga] 


a[AU] 


M p [M Jup ] 


R P [Rjup\ 


Pp IPjup] 


£ = 10° 


£ = 45° 


£ = 90° 


HD 209458 b 1 - 2 


4-7 


0.047(3)** 


0.657(6) 


1.320(25) 


0.28 


1.8 x 10 26 


3.9 x 10 27 


1.2 x 10 28 


OGLE-TR-56 b 3 ' 4 ' 5 


2-4 


0.0225(4) 


1.24(13) 


1.25(8) 


0.63 


4.1 x 10 28 


9.0 x 10 29 


4.0 x 10 30 


HD 189733 b 6 - 7 


0.5-10 


0.0313(4) 


1.15(4) 


1.154(32) 


0.75 


9.6 x 10 26 


2.1 x 10 28 


6.3 x 10 28 


XO-lb 8 


2-6 


0.0488(5) 


0.90(7) 


1.30(11) 


0.41 


1.3 x 10 26 


2.3 x 10 27 


6.7 x 10 27 



** numbers in parenthesis represent the uncertainties on the last digits. *age of the parent star. References: [1] 
Winn et al.(2005); [2] Knutson et al.(2006); a lower R p = (1.26 ± 0.0&)R Jup radius has been recently esti- 
mated by Richardson et al.(2006).[3] Torres et al.(2004); [4] Bouchy et al.(2005a); [5] Santos et al.(2006); [6] 
Bouchy et al. (2005b); [7] Bakos et al. (2006); [8] McCullough et al. (2006); a lower R p = (1.184±0.020)i? y „ p 
radius has been found by Holman et al.(2006). The most recently discovered transiting and potentially inflated 
"hot Jupiters", Tres-2, Hat-P-1, and WASP-1, have not been taken into account. For tidal heating estimates, 
the eccentricity is zero. The dissipation factor <2/£ 2 is set to 10 6 but may vary by several orders of magnitude. 

orbit or non-zero obliquity was also investigated. Since tides circularize the orbit and affect the 
obliquity on timescales (respectively ~ 10 9 yr and ~ 10 6 yr) that are shorter than the system age, 
a continuous tidal dissipation requires a mechanism to maintain the eccentricity or obliquity in 
a non-zero value. Bodenheimer et al. (2001, 2003) calculated that a small eccentricity e 0.03 
forced by a companion would provide the power needed to inflate the HD 209458b's radius, but 
detailed observations suggest that current eccentricities of hot Jupiters are probably smaller (e.g. 
Laughlin et al. 2005). Recently, Winn & Holman (2005) suggested that hot Jupiters may have been 
trapped and maintained in a Cassini state with a large obliquity since the early tidal despinning. An 
accurate estimate of tidal heating is then required to evaluate the role of tidal dissipation inside "hot 
Jupiters". Moreover, it is important to test whether these scenarios are realistic or not. However, 
all previous estimates of tidal heating have considered that the rotation period of "hot Jupiters" 
is synchronous with their orbital mean period. Such an hypothesis is widely justified for most of 
the tidally-evolved solid satellites (e.g. Peale 1999; Wisdom 2004), but we show that this is quite 
improbable for tidally-evolved gaseous planets with a non-zero eccentricity or obliquity. We esti- 
mate the rate of tidal heating within such a non-synchronous "hot Jupiter" and apply it to "inflated" 
planets. We also examine the possibility that the planet's obliquity could have been trapped and 
maintained in an isolated spin-orbit resonance. 



2. Tidal heating within "hot Jupiters" 



We consider the gravitational tides raised by the host star on a short-period planet and follow the 
traditional "viscous" approach of the equilibrium tide theory (Darwin 1908). The secular evolution 
of the spin is then given by (e.g. Hut 1981; Neron de Surgy & Laskar 1997) 

doj K 



dt 



Cn 



(l +x 2 )ci(e)- -2xN(e) 



(1) 
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with 

1 + 3e 2 + f e 4 1 + ^e 2 + %e 4 + ±e 6 

Q(e)= and N ^= (1-^)6 - 

where w is the planet's rotation rate, e its obliquity (the angle between the equatorial and orbital 
planes), x = cos e, <? the orbital eccentricity, and n the orbital mean motion. C is the polar moment 
of inertia and 



K= 3 -^ 
2Q„ 



(GM 



2\ 



MA (R 



M„ 



6 



n, (3) 



where k-i is the potential Love number of degree 2, M* the star's mass, M p the planet's mass, R p its 

radius, and a the semi-major axis. The annual tidal quality factor is Q n = (n Af) where Af is the 

constant time lag. If a separate mechanism maintains a non-zero eccentricity or obliquity, the spin 

evolves to its equilibrium rotation rate (setting dto/dt = 0) 

N(e) 2x 
Q(e) 1 + X 2 

with a time scale r = Cn 2 /K close to ~ 10 5 yr for typical HD 209458b's parameters and 2-values 
comparable to that inferred for Jupiter, ~ 6.10 4 < Q < 2.10 6 (Yoder & Peale 1981). It ends with a 
synchronous rotation if eccentricity and obliquity are simultaneously zero. If not, this equilibrium 
can be higher or lower than the synchronous rotation rate. For solid bodies (satellites, terrestrial 
planets) that have a non-axisymmetric shape, capture in a resonant state where the ratio between co 
and n is commensurable, is also a possible outcome of the spin-down evolution. When eccentricity 
is small, synchronization is the most probable final state (e.g. Goldreich & Peale 1966), but other 
non-synchronous resonances are possible with a large orbital eccentricity (e.g. Correia & Laskar 
2004). The capture probability is generally an increasing function of the asymmetry parameter 
(B - A)/C, where A and B are the equatorial moments of inertia. However, determination of the 
second-degree harmonics of Jupiter's and Saturn's gravity field from Pioneer and Voyager tracking 
data (Campbell & Anderson 1989; Jacobson 2002) provided a crude estimate of the (B - A)/C 
value that is lower than ~ 10~ 5 for Saturn and ~ 10" 7 for Jupiter. These values are more than one 
order of magnitude lower than the Moon's or Mercury's value, leading to insignificant chances of 
capture. Furthermore, the detection of this asymmetry is questionable. If such an equatorial bulge 
originates from local mass inhomogeneities driven by convection, it is probably not permanent and 
must have a more negligible effect if averaged spatially and temporally. In this context, we assume 
that the tidally-evolved rotation state of "hot Jupiters" with a non-zero obliquity or eccentricity is 
given by Eq.©. Once this equilibrium is reached, tidal energy continues to be dissipated in the 
planet at the expense of the orbital energy so that EtUe = _ (GM p M*)/ (2 a 2 ) xda/dt with (e.g. Hut 
1981; Neron de Surgy & Laskar 1997) 



da o I K 
— = 4-a 



dt XGM+Mp 



N(e)x--N a (e) 
n 



(5) 



where N a (e) = (1 +31e 2 /2 + 255e 4 /8 + 185e 6 /16+25e 8 /64)/(l - e 2 ) 15 ' 2 . Thus, with the equilibrium 
rotation given by Eq.(@), we obtain: 

N 2 {e) 2x 2 



Elide — IK 



N a (e) ■ 



(6) 



Q(e) 1 +x 2 
At second order in eccentricity, we get 
2 K 

E tide (e, e) = 5- [sin 2 e + e 2 (7 + 16 sin 2 e)l , (7) 

1 + cos 7 e L v n 
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Fig. 1. Rate of tidal dissipation within HD 209458b as a function of eccentricity for a 45° (solid 
thin line) and 90° (solid thick line) obliquity. The synchronous case is plotted with dashed lines for 
comparison. The dissipation factor Q„/k2 is set to 10 6 . 

which is always larger than in the synchronous case E s t j^'(e,s) = K (sin 2 e + le 2 ) (e.g. Wisdom 
2004). In FigQ] we compare the rate of tidal heating within a non-synchronous and synchronous 
planet as a function of the eccentricity (0 < e < 0.25) for two different obliquities. The ratio 
between them is an increasing function of both eccentricity and obliquity. For e ^ 0, as observed 
for "hot Jupiters", it reaches ~ 1.5 and 2 at 45° and 90° obliquity respectively, slightly increasing 
at larger eccentricity. In Table. 1, we compute the rate of tidal heating within each inflated "hot 
Jupiter" at 10°, 45°, and 90° obliquity. It reaches ~ 10 29 -10 3() erg.s 1 for OGLE-TR-56b and HD 
189733b due to their proximity to the central star. Using the evolutionary models of Chabrier et 
al.(2004), we estimated that an additional power of 0.5 - 1 .5 x 10 27 erg.s -1 is required to explain the 
present inflated radius of each planet (if coreless). Table. 1 shows that such values can be provided 
by obliquity tides. 

3. Capture in isolated spin-orbit resonances 

Winn & Holman (2005) suggest that obliquity tides can be maintained if the planet has been trapped 
in a spin-orbit resonance since the early despinning process. Their study uses the traditional defini- 
tion of Cassini states, which requires that the nodal precession frequency and the inclination remain 
constant in time. Here, this hypothesis is relaxed and we consider the effects of gravitational pertur- 
bations caused by other fellow planets or a protoplanetary disk, which induce secular variations in 
orbital parameters of a "hot Jupiter". In this context, the precession equations for the obliquity and 
the general precession tfr can be written as the sum of quasi-periodic terms (e.g. Laskar & Robutel 
1993; Correia & Laskar 2003): 




(8) 



k 




(9) 



k 
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a=0.5AU a=0.1 AU a=0.05AU 




0.1 1 10 0.1 1 10 0.1 1 10 

Rotation period (days) Rotation period (days) Rotation period (days) 

Fig. 2. Obliquity capture probabilities in resonance (v* = -10.0 "/yr) under the effect of gravita- 
tional tides, as a function of the rotation period at a) 0.5 AU b) 0. 1 AU c) 0.05 AU. Jk is the ampli- 
tude of secular orbital perturbations (see. EqJSJl. Physical parameters are taken from HD 209458b. 



where are secular frequencies of the orbital motion, Jk and cpk the respective amplitude and phase, 
a — (3 n 2 Ed) I {2 co) is the precession constant, and Ed = (C —A)/C is the dynamical ellipticity. For 
a fluid planet in hydrostatic equilibrium, Ed = (k / Rp)(3GC) x to 2 = Ko u> 2 where kf is the fluid 
Love number (~ 0.93 for Jupiter), so that 

a = (3/2)K n 2 oj oc cj/a\ (10) 

In the following, we do not consider the variations in frequencies caused by the evolution of 
the protoplanetary disk or planetary migrations, and do not explain the details of the origin and 
calculation of these frequencies. For small amplitude variations of the eccentricity and inclination, 
typical values of Jk are low (\Jk\ <s |vjt|), and the condition of resonance between the spin axis 
precession and the secular frequency v* yields 

a cos e sine + Vk sine Jk cose , (11) 

which can be seen as a generalized condition for Cassini states and which has 2 or 4 roots for each 
isolated frequency Vk- For low values of \Jk\, the stable states can be approximated well by 

Jk 

tanei ± ; cose2 — —Vk/a . (12) 

Vk + a 

Unless a = |vj;|, the state 1+ is close to 0°, therefore producing negligible heating and state 1_ is 
close to 180° and thus unreachable by tides. We then focus only on state 2, which can maintain a 
significant obliquity. To test the possibility of capture in such a state, we consider a simple scenario 
where a "hot Jupiter" forms at a large orbital distance (~ several AUs) and migrates inward to 
its current position (~ 0.05 AU). Before the planet reaches typically ~ 0.5 AU, tidal effects do not 
affect the spin evolution, but the reduction in the semi-major axis increases the precession constant, 
so that the spin precession may become resonant with some secular orbital frequencies Vk - This kind 
of passage through resonance generally causes the obliquity to change (e.g. Ward 1975; Laskar et 
al. 2004) raising the possibility that the obliquity has a somewhat arbitrary value when it attains 
~ 0.5 AU. At this stage, tidal effects can increase the obliquity toward the evolving equilibrium 
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value computed from Eq.© (setting ds/dt = 0) 

In Me) 2n , ,\ 
cosfe,) = —±{ - — 1 + 6e 2 ) , (13) 

which tends to 90° for initially fast rotation rates {to » «). As the rotation rate is decreased by tides, 
the equilibrium obliquity is progressively reduced to zero degrees. It is then possible that the obliq- 
uity crosses several secular spin-orbit resonances (one for each frequency Vk) in both ways (increas- 
ing and decreasing obliquity) and that a capture in one of those occurs. Inside the resonance island, 



the restoring torque causes the obliquity to librate with amplitude Ax = A cos s — 4 V-Zt sin £2/0: 
around £2 (e-g- Correia & Laskar 2003). The probability of capture in this resonance can be esti- 
mated using the analytical approach of Goldreich & Peale (1966). When a linear approximation of 
the tidal torque (Eqf2]i is used around the resonant obliquity, we have 



^ca P = r and T) = 

1 + 77 1 



(1 - 3 cos £2) - + 2 cos £2 



Ax . (14) 



n sin 2 £2 (2 - cos £2 a>/ri) 

In Figf2] we have plotted the capture probabilities at 0.05, 0.1 and 0.5 AU as a function of the 
rotation period for the different amplitudes (J^) and frequencies (v*) characteristic of the Solar 
System (e.g. Laskar & Robutel 1993). As a reasonable example, we choose Vk = — 10"/yr, but 
our results are not affected by changes in this value. Assuming an initial rotation period of ~0.5 
day, the capture is possible and even unavoidable if Jk > 0.1"/yr at 0.5 AU. On the contrary, 
the chances of capture at 0.05 AU are negligible (< 1%) because a decrease in the semi-major 
axis leads to an increase in the precession constant and reduces the width of the resonance. We 
compared these theoretical estimates with numerical simulations. For that purpose, the obliquity 
equations were integrated in the presence of tidal effects (Eqsf2][8]and|9]) considering 1000 initial 
precession angles, for each initial obliquity, equally distributed over - 27r. Statistics of capture 
were found to be in good agreement with previous estimates. Once captured, both spin-down and 
migration processes continue to slowly modify the precession constant value until the equilibrium 
rotation rate (Eqj4j) and semi-major axis are reached, respectively, on typical 10 5 -10 6 yr and 10 5 - 
10 7 time scales. To test their influence on the capture stability, we again performed some numerical 
simulations for various initial obliquities and secular perturbations over typically ~ 5.10 7 yr. The 
migration process was simulated by an exponential decrease in the semi-major axis towards 0.05 
AU with a 5.10 s - 10 7 yr time scale. The obliquity librations were found to be significantly shorter 
than spin-down and migration time scales so that the spin trajectory follows an "adiabatic invariant" 
in the phase space. Then, the average obliquity tracks the evolving value cos £2 = -v^/a without 
leaving the resonance. Under these conditions, using Eqs.fHJi, ( fTUl ), and ( TTZl i. the theoretical final 
obliquity is then given by 



COS (S final) = A / y— , (15) 

where 11 = (-vt)/(3 K^n 3 ), which does not depend on the rotationrate. At 0.05 AU, we have u <sc 1 
so that Sfi na i =a 90° and the final spin frequency is <±> eq 2 -\furt. However, Eq.(f2| indicates that 
the tidal torque dramatically increases both with spin-down and inward migration processes (s oc 
a -i5/2 w -i-) jf t n j s tidal torque exceeds the maximum possible restoring torque (Eq[8]l, the resonant 
equilibrium is destroyed (the evolution is no longer adiabatic). For a given semi-major axis, the 
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stability condition then requires that the rotation rate must always be higher than a threshold value 
o)crit.t which is always verified if co eq > (jj cr j t .. We found that this condition can be simply written as 

tan(e f ina i) < J k x r (16) 

where t = C n 2 /K is the time scale of tidal despinning. It then follows that the final obliquity of 
the planet cannot be too large, otherwise the final equilibrium rotation rate would be too low (EqHJi 
and the obliquity quit the resonance. For instance, taking = 1 "/yr at 0.05 AU (our highest value 
in Figf2]i, we need e/, m / < 21°. This value drops to e/,, M ; < 2° for the more realistic Jk = 0.1 
"/yr amplitude. However, such a low resonant obliquity at 0.05 AU is highly unlikely because, 
according to Eq.dTSt, it requires very high values of the orbital secular frequencies |vt| > 7.2 x 10 6 
"/yr. At 0.5 and 0.1 AU, critical obliquity values are respectively ~ 90° and 83.4° for Jk — 1 "/yr, 
89.9 and 41° (for Jk - 0.1 "/yr) and require more reasonable orbital secular frequencies so that a 
stable capture is possible. We empirically retrieve the stability criteria ( fToT ) for the final obliquity 
with excellent agreement in our simulations. When the obliquity leaves the resonance e%, it rapidly 
decreases towards the stable state e\ + =t 0°, producing negligible obliquity tides (EqfT2b. 

4. Conclusions 

We computed tidal heating within non-synchronous "hot Jupiters" and found that obliquity tides 
may quantitatively explain their anomalously inflated radius. However, locking a "hot Jupiter" in 
an isolated spin-orbit resonance seems highly unlikely at 0.05 UA, although this could be possible 
for giant planets at larger distances (typically > 0.1 AU). We did not consider the overlapping of 
several spin-orbit resonances causing the obliquity motion to become potentially chaotic (Laskar & 
Robutel 1993). A more realistic description of perturbations in a protoplanetary disk or extra-solar 
systems is needed to examine these scenarios. 
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